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Stark effect upon the effective mass and radius 
in a tight-binding exciton model. 
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With a Green's function formalism we obtain the eigenvalue spectrum of a tight-binding one- 
dimensional exciton model characterized by a contact interaction, a Coulombic electron and hole 
attraction, the Heller- Marcus exciton-hopping energy and an external constant and homogeneous 
electric field. The resulting eigenvalue spectrum, in the form of an unevenly spaced Wannier-Stark 
ladder with envelope profiles, is used to obtain the effective mass of the exciton by the application 
of the Mattis-Gallinar effective mass formula [D. C. Mattis and J. -P. Gallinar, Phys. Rev. Lett. 53, 
1391 (1984)]. We obtain positive and negative effective masses for the exciton. The inverse effective 
mass may oscillate periodically as a function of the inverse of the electric field, with the frequency of 
oscillation linearly dependent upon the tight- binding hopping matrix element. The exciton radius is 
also obtained with the Green's function formalism, and it too exhibits Keldysh-like field dependent 
oscillations, as well as abrupt variations associated to strongly avoided crossings in the eigenvalue 
spectrum. Finally, some comments are made about the experimental relevance of our results. 
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I. INTRODUCTION 



As shown by Mattis and Gallinar |y, 0, the mass of 
an exciton is not, in general, given by the intuitive result 
of the sum of its electron and hole masses, but it may 
depend upon excitonic internal quantum numbers, such 
as, for example, the internal kinetic energy of the exci- 
ton 0| , or the Heller-Marcus exciton-hopping energy || . 
Beautiful experimental confirmation of some of these sur- 
prising, and interesting, predictions, was given by Cafolla 
et al. |Q in 1985 from electron-energy-loss spectroscopy. 
Up to now, however, no detailed study or analysis of a 
specific exciton model || has been made p8| in the light 
of the Mattis-Gallinar formula (MGF) (cf. Eq.f])) and 
of its interesting counterintuitive predictions. In view of 
this, and to further study these predictions, here we shall 
present a specific, simple, one-dimensional exciton model 
susceptible to an exact analysis for its eigenvalue spec- 
trum, even while in the presence of an external electric 
field. Thus enabling the study (in an exact manner) of 
diverse aspects related not only to the mass of an exciton 
through the application of the (MGF), but also of those 
related, for example, to another important general exci- 
ton property, such as its radius. As we shall see, some of 
the shared characteristics of these two exciton parameters 
(mass and radius) can be understood in terms of the un- 
derlying eigenvalue spectrum of the exciton, which, due 
to the applied external electric field, is obtained to be 
in the form of an unevenly spaced Wannier-Stark ladder. 
This ladder, in turn, exhibits Keldysh-like oscillations || 
and strongly avoided crossings of energy levels, with dra- 
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matic consequences for the mass and radius of the exciton 
worthy of close examination. 

The plan of this paper is then as follows: in Section || 
we present the mathematical Green's function formalism 
to be used to obtain the eigenvalue spectrum, as well as 
the different parameters that go into the exciton model 
to be considered. In Section III we exhibit a collection of 



graphical results, and in Section IV conclusions are pre- 
sented and discussed. Finally, an Appendix with several 
analytical results is included. 



II. EXCITON MODEL 

The one-dimensional exciton model that we will con- 
sider || , akin to that of the Merrifield (7) exciton, can be 
defined in terms of the matrix elements of the excitonic 
Hamiltonian TL. If |n, m) represents a Wannier orthonor- 
mal basis states in which the electron is localized at site 
n, and the hole at site m, the nonvanishing matrix ele- 
ments of H will be taken to be 



(n, m\H\n, m) = W + V(n — m), 



(1) 



where W represents the combined half-bandwidth of the 
electron plus hole, n — m is the relative electron-hole co- 
ordinate and V(n—m) is the effective two-body potential 
through which the electron and hole interact with each 
other Q], and with the external electric field, if present 
p9[ . We will assume that V(n — m) is given by 



V(n) = -V 8 n , c 



(1 - <W>) 7~~T 
n 



(2) 
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where n = 0, ±1, ±2, V D measures the contact interac- 
tion j7j, or electron and hole attraction when they are 
at the same site on top of each other, while, on the 
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other hand, V measures the long-range part of the 
attraction, modelled here by a Coulombic l/\n\ tail. Fi- 
nally, the interaction of the exciton electric dipole mo- 
ment with the external electric field E is given by an, 
with a — ae E, where "a" is the lattice constant and 
(— e) the electron charge. 

The non-diagonal matrix elements of TL associated to 
the kinetic energy of the exciton are given by 



(n, m\TL\n, to ± 1) = — 1%, and 



(n, m\Tl\n ± 1, to) 



-U 



(3) 



With th representing the usual valence-band width pa- 
rameter associated with nearest-neighbor jumping of the 
hole, and t e the corresponding conduction-band width 
parameter associated to the electron. Thus, W = 2(th + 
t e ). Finally, we introduce the exciton- hopping (or Heller- 
Marcus) matrix elements given by 



(m, m\TL\n, n) — H(n — to), n =/= m 



(4) 



whereby the entire (Frenkel-like) exciton jumps as a 
whole entity between sites n and m, with the amplitude 
H(n — m). 

According to the (MGF) jl], ||, the translational mass 
M„ of an exciton can be written as 



m, 



m h 



1 



Kn 
W 



m h H n 



(5) 



H, 



where m e and nrih are, respectively, the electron and hole 
effective masses, K n and H n are the expectation values 
of the kinetic and exciton-hopping energies of the exci- 
ton in the nth bound or localized state of the interaction 
potential, and, finally, Mp and Hp are the values respec- 
tively taken by M n and H n for a Frenkel-like, strongly 
localized exciton. 

As Eq (|) shows, the mass M n of the exciton can 
be determined once the explicit values of K n and H n 
are given. These values can be obtained through the 
Hellmann-Feynman theorem, by writing || 



K n = 



dE n 
dX 



, and 



dE n \ 



(6) 



where E n = E n (k — 0) is the eigenvalue of TL in the nth 
bound or localized state of the exciton at k = 0, with k 
being the center-of-mass wave vector of the exciton jj). 
In Eq. (||), E n is made to depend upon the dimensionless 
parameters A and fj, through the following simple scaling 



th — > Xth 



and 



H(n — to) — > j.iH{n — to). 



(7) 



A Green's function formalism j|, permits then the 
evaluation Q of the energy eigenvalues E n , by finding 
the simple poles z n (m) of the continued fraction G m ^ m (z) 
given by 

e 2 

Gj m {z) = z-V(m) 



z - V(m + 1) 



z — V(m — 1) 



(8) 



Where e — A (t e +th), and to = 0, ±1,±2, ... is a fixed 
integer. Due to the peculiar nature of the Heller-Marcus 
pOf exciton-hopping energy ||, the Term V(0) appear- 
ing in Eq. (|) is shifted so that § V(0) -> V(0) + fiH F . 
The set of all eigenvalues E n of TL can then be found and 
ordered appropriately || from the relationships E n — 
z n (m) + W , in a manner over-all independent of the cho- 
sen m's (from now on for simplicity we write z n (m) = z n ). 

The mean-square radius of the exciton can also be 
found in terms of G~ 1 m (z). In fact, the mean-square 
radius of the exciton is defined in terms of the relative 
electron-hole coordinate to, as 



Rl= J2 CL 2 m 2 \F n (m)\ 2 , 



(9) 



where F n (to) is the normalized relative coordinate wave- 
function of the exciton || . As shown in the Appendix, 
the Green's function formalism Q permits then the di- 
rect evaluation of |F„(to)| 2 , obtaining finally 



m— +0O 

R n = 

m— — oo 



(10) 



In contrast with the evaluation of the set of eigenvalues 
E n , all values of to are seen in Eq ([l0]) to be in principle 
relevant for the evaluation of i? 2 . In the next section 
3 we shall display some of the results obtained through 
the use of the above described formalism applied to our 
exciton model. 



III. RESULTS 

In this section we display three types of results: en- 
ergy spectra, effective masses and root-mean-square radii 
of the exciton. Figures (1) and (2) were obtained start- 
ing from the numerical evaluation |L0| of the continued 
fraction in Eq.(||). In Fig (1), to was chosen for eval- 
uation purposes as to = due to absence of Coulomb 
interaction, while in Fig (2) m's different from zero were 
appropriately chosen for ease of computation. Because of 
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the contact interaction V being V /a = 10, level n = in 
the ladder in Fig (1) is moved down to the level n = — 10 
where a degeneracy occurs at t = 0. While in Fig (2), 
with Coulomb interaction, no such degeneracy occurs for 
the selected parameters. In both Figures (1) and (2), sev- 
eral (strongly) avoided energy crossings are evident and 
are associated to the envelope profiles. In Fig (1), the 
envelope profile which starts at z n /a — —10 is that of 
the single bound state |ll| of the contact interaction V 
(in absence of the electric field), i.e., 

E n - W = -\Jv o 2 + W 2 . (11) 



Inhomogeneous Wannier-Stark ladder 

Vo/a = 10, V/a=0 

3 T 




2t/a 



Inhomogeneous Wannier-Stark ladder 

Vo/V=3, Vo/a=10 



6 T 




2t/a 



FIG. 2: Inhomogeneous Wannier-Stark Ladder energy levels 
E n — W for the exciton in terms of z„(m)/a vs 2t/a, with 
V /V = 3 and V a /a = 10; n runs from —10 to 6, the chosen 
m's are not all zero and t = t e = th- Note the increment in 
the degree of flatness down in the ladder where a large num- 
ber of derivatives is vanishingly small at t = 0. Noteworthy 
characteristics of the ladder are: oscillations as t increases, 
strongly avoided energy crossings and envelope profiles. 



FIG. 1: Inhomogeneous Wannier-Stark Ladder energy levels 
E n — W for the exciton in terms of z n (m)/a vs 2t/a, with 
V = and Vo/a = 10; n runs from —13 to 2 (skipping n = —1 
for visual convenience), the chosen m is zero and henceforth 
t = t e = th- Note the increment in the degree of flatness down 
in the ladder where a large number of derivatives is identically 
zero at t = 0. Noteworthy characteristics of the ladder are: 
oscillations as t increases, strongly avoided energy crossings 
and envelope profiles. Energy levels n — and n — — 10 are 
degenerate at t — 0. 

As for Fig (3), it was numerically obtained using a 
scaling property Q (valid when Hp = 0) related to the 
energy E n (k;t), namely, 

E n (k;t) - W = z n (k;t) = z„(0; t cos ka). (12) 

Different exciton masses are to be associated to levels 
n = and n = —10, which have respective curvatures of 
different magnitude at k = in Fig (3). 



Figures (4), (5) and (6) display a small representative 
sample of our results for the effective exciton mass M n 
p"2[ . Interestingly enough, we find both positive and neg- 
ative effective masses with | m "^ mh \ < 1 when Mp = oo 
(i.e., in absence of the Heller-Marcus effect in Figures 
(4) and (5) ). In Fig 6 we find that the exciton-hopping 
effect U shifts the mass away from negative values but 
still allows for effective negative masses. That the effec- 
tive mass of the exciton may indeed become negative in 
the presence of an external electric field is an interest- 
ing result found here that has not been properly noted 
fl3f in the previous literature ||, [|, ||, [|, and might be 
of experimental relevance for the excitonic Stark effect. 
The noticeable visual envelope for the inverse mass in Fig 
(4) arises from the bound state of Eq.(0), or envelope 
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Wannier-Stark ladder levels n=0,-10 

Vo/a=TO, VorV=3,t/a=5 
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FIG. 3: Plot of Wannier-Stark ladder levels n = (red) and 
n = — 10 (blue) vs ka in the first three Brillouin Zones. The 
strongly avoided crossing of corresponding levels shown in fig- 
ure |^ appears here close to the borders of the Brillouin Zones. 
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FIG. 4: Plot of the inverse effective dimensionless mass (m e + 
rrih)/M n vs 2t/a, with V /a — 10 and Mf — oo, for the 
energy levels n — 0,-11,-12,-13. Both here and in the 
Wannier-Stark ladder in figure these levels give rise to a 
noticeable visual envelope and abrupt changes in the exciton 
mass. The amplitude of the oscillations decreases for large 
values of 2t/a. 



profile of Fig (1). It is given by [O 



m h 



M n 



vT 



I''- 



(13) 



where x = W/V a . 

In Figures (5) and (6), the inverse of the mass is seen 
to oscillate (for cT 1 — > oo) as a function of 1/a. The 
period P of these Keldysh-like oscillations can be shown 
to be (see Appendix) 



4 V t 



(14) 



Numerical results |Q furthermore show Eq. ( |l4| ) to be 
still valid when the Coulomb interaction V is present, 
i.e., the period of these oscillations is independent of V. 

Our three remaining figures display typical behavior 
of the root-mean-square radius R n of the exciton in our 
model. Figures (7) and (8) display an oscillatory behav- 
ior as R n diverges with increasing t — > oo. Superimposed 
on both figures (7) and (8) is an overall hyperbolic "pat- 
tern" , given by the formula (see Appendix) 



Rn 

a 



(15) 



which expresses the value of R n under the action of the 
electric field alone, i.e., for a homogenous equally spaced 
Wannier-Stark ladder, with eigenvalues E n — W = na. 

In Fig. (8) the radius R n /a for n = starts at zero, and 
for small enough 2t/a it increases linearly (cf. Eq.djlj)) 
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FIG. 5: Plot of the inverse effective dimensionless mass (m e + 
m h )/M n vs Vo/a, with t/V = 0.7, V/V = and M F = oo, 
for two energy levels (n = —1,-2) from the Wannier-Stark 
Ladder in figure |l| As a -1 — > oo, the curves become periodic 
in the inverse electric field 1/a, with period P given by P = 
j{V /t) w 1.122 . 



until it reaches the value of 2t/a ~ 1. At this point 
one observes an abrupt variation of R n /ce for both levels 
(n=0, -10), and a simultaneous strongly avoided crossing 
in the Wannier-Stark ladder. From then on the contri- 
bution to the linear increase in the radius (cf. Fig. (2)) 
is given by the following level down in the Wannier- 
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Inverse of the Exciton mass 
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FIG. 6: Plot of the inverse effective dimensionless mass (m e + 
rrih)/M n vs V /a for the energy level n = —2, with t/V = 
0.2, V/Vo = 0, and three values of the Wannier-Frenkel ratio 
(m e + m,h) I 'Mf =wf. The Heller Marcus effect shifts the mass 
away from negative values. As a -1 — » oo, the curves become 
periodic in the inverse electric field 1/a, with identical period 
P given by P = f (V D - H F )/t » 3.927 . 
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FIG. 8: Plot of the dimensionless exciton radius R n /a vs 2t/a 
for two levels n = (red) and n = —10 (blue), with V /a = 10 
and V /V = 3. The exciton radius shows abrupt variations 
due to the strongly avoided crossings in the corresponding 
energy spectrum (figure ^ . Note the envelope associated to 
each one of the exciton radii, arising from the envelope in the 
Wannier-Stark ladder. 
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FIG. 7: Plot of the dimensionless exciton radius R„/a vs 
2t/a for three levels (n = —1,-3,-4), with V /ot = 10 and 
V /V = 3. At t = exciton motion is frozen, with expected 
inter-particle distance equal to \n\. As t — > oo exciton radius 
diverges in an oscillatory manner. 



Stark ladder (here n=-10). These abrupt variations in 
the exciton's size, exemplified here with these two lev- 
els (n — 0,-10), arc, however, a general feature of our 
model jl^]. One which will always occur associated to 
a strongly avoided crossing in the respective eigenvalue 
spectrum, and one which has no counterpart in the study 
of, for example, the standard Stark effect for the hydro- 



FIG. 9: Plot of the dimensionless exciton radii R n /a for n — 
(red) and n = — 10 (blue) vs ka in the first Brillouin Zone. 
The strongly avoided crossing of corresponding levels shown 
in figure appears here as abrupt variations of the radii close 
to the borders of the first Brillouin Zone. 

gen atom |6| . 

Finally, in Fig. (9) we display the fc-dependence of 
Rn/a for the same levels as in Fig. (8) at t/a = 5. The 
abrupt variations in Fig. (8) are displayed in Fig. (9) 
close to the first Brillouin Zone borders, and oscillations 
of R n /a in one curve are transposed onto the other. 

Having displayed graphically some of our extensive [T^ ] 
numerical results, in the next section some conclusions 
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are presented and comments made about the relevance 
of our exciton model. 



IV. CONCLUSIONS 

Although simple, the exciton model considered previ- 
ously displays quite interesting behavior. Most notably, 
oscillating (negative) excitonic masses and abrupt varia- 
tions in the exciton's size (among others). To the best 
of our knowledge, these are novel predictions that have 
not been considered before in the literature, nor have 
these effects been experimentally observed. Since the 
strongly ionizing influence of an external electric field 
upon a weakly bound exciton is well-known Jl3| ], it is 
clear that the difficulty of experimental observation of 
some of these effects might be a large one. Furthermore, 
it is unclear to what extent the previous theoretical re- 
sults may be modified by, for example, the inclusion of 
a more realistic exciton (band) dispersion law, as well as 
the lifting of the one-dimensionality of the model. To 
this extent, further investigation of our model's interest- 
ing predictions is warranted in the future. 

From a self-consistent point of view it is worthwhile to 
point out that, within the context of the model's given 
tight-binding dispersion, the semiclassical effective po- 
tential |i~4| V f.jf.(n) associated to V(n) (cf. Eq. (^)), is 
given by^l 



v *ff. ( n ) = 



V 2 (n) — a 2 EV(n) 



(16) 



(where E is the semiclassical energy of the exciton and 
fi* is the reduced mass of the electron and hole). Since 
V e ff. (n) is positive and diverges for |n| — » oo, this pre- 
cludes the occurrence of the well-known Jf3[ tunneling 
through the potential barrier that appears in the stan- 
dard Stark effect for the hydrogen atom, or for a non 
tight-binding exciton in the parabolic effective mass ap- 
proximation [ fi"3| . Thus, our bound exciton states on the 
lattice are truly localized Jl6|], instead of being simply 
metastable ones as for a hydrogen-like atom on the con- 
tinuous manifold and consequently, both the mass 
and the radius of the exciton are well-defined concepts 
for these states. It is to be noticed, finally, that although 
presented in the context of the (MGF), our exact results 
for the effective mass M„. of the exciton have also been 
obtained independently Q from the exciton dispersion 
law E n (k). Thus the validity of the (MGF) has in this 
respect also been checked. 

In summary, we have considered a very simple one- 
dimensional exciton model in the presence of a static and 
homogenous electric field, and this model has been found 
to display highly interesting behavior for some of its prop- 
erties, such as mass and radius Further theoretical 
work should then be able to elucidate to what extent 
some of these predicted effects might be experimentally 
observable. 



APPENDIX 

In this Appendix we derive several analytical results 
used in this paper. First, we obtain Eq.(flO|) in the text 
from the definition in Eq.(^). Based on the Green's func- 
tion formalism 0, the continued fraction G m . m (z) can 
be written as: 



Gr, 



(A.I) 



where the z m ''s ar e sim ple poles of G m ^ m (z). Multiplying 
both sides of Eq. ( |A.f| ) by (z — z n ), we obtain: 



(z — z n ) G m ^ m (z) — S ■ 



^\F m ,{m)\\ (A.2) 

and by taking the limit (z — > z n ), we get 

lim (z - z n ) G m>m (z) = \F n (m)\ 2 , (A.3) 

z — *z n 

or equivalently 

[z - z n ) _ 1 



lim 



z^Zn Grri,m(z) -q£ {Gm,m{z)) , 



= \F n {m)\\ (A.4) 



By substituting the result of Eq.(A.4) into Eq. (J9|) we ob- 
tain Eq.©. 

Now we obtain Eq. (|lj) which gives the period P of the 
Keldysh-like oscillations observed in the energy and the 
inverse mass of the exciton in the limit t/a — ► oo. On one 
hand, we start with Eq.(||) written in the absence of the 
Coulomb interaction, i.e., with V(m) = — V S mt o + a m 
p3, and then we use the following relationship 



(I) 1 



y Jv-i{y) 

2v J u {y) 



- 1 



(I)' 



(A.5) 



(i/ + l)(i/ + 2) 



{v + 2){v + Z) 



which expresses Bessel functions J v {y) of the first kind 
in a continued fraction form | oL Using the following 
identities 



Ju+i{y)J-u{y) + Ju{y)J- v -i(y) 



-2 sin(i/7r) 



Jv-i{y) + Jv+\{y) = — J v (y), and (A.6) 

y 

J v {-y) = {-l)"J v {y), (A.7) 
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Eq. (g) can then be written (for m = and e = 2t) as 

V , (A.S 



x a sen(yir) 

u o,oK z ) - —~ 



it J v (y)J^ v (y) 

where v = z/a and y = At/ a. 

The eigenvalues E n are obtained from the following 
implicit equation upon the z n 's, namely 



sin 



(^) 



V ir 



J -Mr, 



obtained in turn from G (z) — 0. When a 
use the asymptotic Bessel functions form jL7) 

lim J± v {y) = 



y^ca 



2/1 1 

COS U T — ^7T 7T 

Try V 2 4 



(A.9) 



co, we 



(A.10) 



Upon substitution of Eq.( A.ld| ) into Eq.(A.G), we obtain 
after some algebraic manipulations that 



it . (z n 
— sin 7T 

Vo U 



> / 1 z n 
COS 7T 

2 a 



1 I 1 Z « 
sin 7T 

2 a 



,'4A (At 
2 sin I — cos I — 
a ) V a 



— sin — 7r = cos — 7r + sin — — : A. 11 
A V a J \ a / \A a J 

z n /a as a function of V /a will then be periodic with 
period P — jA, where A = ^f-, and so will the inverse 
mass of the exciton according to Eq.(||). 



To end this Appendix we indicate briefly the derivation 
of Eq.@. By setting V a = V = H F = in Eq. (|), we 
find 



7T Vtt / 

° — — m I I J — — +m 1 I 

a - Ui 



(A.12) 



which can be proven |12| in a manner analogous to that of 
Eq. ( |A.S| ). One now evaluates the derivative of G~ 1 m (z) 
appearing in Eq. (|io|), as 



d_ 



G m . m ( z )\z=z n — j (H\ j / 4t V (A. 13) 



where z n — not. Substituting Eq. ( A. 13 ) into Eq. (|1C|), 
there results 



m— +oc f At\ f At 

R 2 n = X! (-l) m+ ™ a 2 m 2 J(„_ m ) f — J J_ (n _ m ) f — 



which reduces to 



2 m— +oo 



m— — oo 



At 



(A.14) 



(A.15) 



by using Eq. (A. 7). Eq. (|15| ) finally ensues after some 
algebraic manipulations ]12[ on Eq. (A.15) 



[1] D. C. Mattis and J.-P. Gallinar, Phys. Rev. Lett. 53, 1391 

(1984) . 

[2] D. C. Mattis, Rev. Mod. Phys. 58, 361 (1986). 

[3] J.-P. Gallinar and D. C. Mattis, Phys. Rev. B 32, 4914 

(1985) . 

[4] A. A. Cafolla, S. E. Schnatterly, and C. Tarrio, Phys. 
Rev. Lett. 55, 2818 (1985). 

[5] J.-P. Gallinar, J. Phys.: Condens. Matter 8, 3103 (1996), 
and references therein. 

[6] J. Callaway, Quantum Theory of the Solid State (Aca- 
demic Press, New York, 1974). 

[7] R. E. Merrifield, J. Chera. Phys. 34, 1835 (1961). 

[8] J.-P. Gallinar, Phys. Lett. 103A, 72 (1984). 

[9] E. N. Economou, Green's functions in Quantum Physics, 
Springer Series in Solid-State Sciences 7 (Springer- 
Verlag, Berlin, 1979). 
[10] W. Press, S. A. Teukolsky, W. T. Vetterling, and B. P. 
Flannery, Numerical Recipes: The Art of Scientific Com- 
puting (Cambridge University Press, New York, 1986). 
[11] J.-P. Gallinar, Phys. Rev. B 36, 1782 (1987). 
[12] J. El-khoury, Master's thesis, Universidad Simon 
Bolivar, Venezuela (1999), available upon request at 
clkhoury@ucla.edu. 



[13] E. I. Rashba, Excitons (North-Holland, Amsterdam, 
1982), edited by E. I. Rashba and M. D. Sturge. 

[14] N. W. Aschcroft and N. D. Mermin, Solid State Physics 
(Holt, Rinehart & Winston, New York, 1976). 

[15] J.-P. Gallinar and J. El-khoury (2001), unpublished. 

[16] S. I. Serdyukova and B. N. Zakhariev, Phys. Rev. A 46, 
58 (1992). 

[17] Handbook of Mathematical Functions (Dover, New York, 
1965), edited by M. Abramowitz and I. A. Stegun. 

[18] With or without an externally applied electric field. 

[19] No external field is considered in M. 

[20] The Heller-Marcus term acts only when the electron and 
hole are at the same lattice site. 

[21] Furthermore, our Green's function matrix elements 
G m ,m{z) have no cuts, i.e., they display only isolated 
poles on the real axis. 

[22] Largely independent of boundary conditions, since it is 
noteworthy that placing our system in a finite "box" (so 
as to make the limit a — » non-singular) does not alter 
significantly any of our physical results |12]. 

[23] When there is exciton-hopping, V is simply replaced by 
V -H F 



